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Abstract. 

We review 'particle partitioning entanglement' for itinerant many-particle systems. 
This is defined as the entanglement between two subsets of particles making up the 
system. We identify generic features and mechanisms of particle entanglement that 
are valid over whole classes of itinerant quantum systems. We formulate the general 
structure of particle entanglement in many-fermion ground states, analogous to the 
'area law' for the more usually studied entanglement between spatial regions. Basic 
properties of particle entanglement are first elucidated by considering relatively simple 
itinerant models. We then review particle-partitioning entanglement in quantum states 
with more intricate physics, such as anyonic models and quantum Hall states. 

PACS numbers: 03.67.Mn, 03.75. Gg, 64.70.Tg, 71.10.-w 
1. Introduction 

Bipartite entanglement in many-particle systems, i.e., the entanglement between one 
part (A) of a system and the rest (B), has grown into a widely studied topic in the 
last few years. Usually, the partitioning is spatial, so that A is a collection of lattice 
sites or is a connected region of space. In this article, we will consider an alternate 
form of partitioning, namely particle partitioning. With the wavefunction expressed in 
first-quantized form, one can meaningfully partition particles rather than space, and 
calculate entanglements between subsets of particles. Since each particle has a label in 
first-quantized wavef unctions, indistinguishability does not preclude well-defined subsets 
of particles. Note that, with such partitioning, A or B do not correspond to connected 
regions of space. Also note that particle partitioning is only defined in itinerant systems 
where the particles hop, and thus has no meaning for pure spin models. 

Particle partitioning entanglement is generally quite different from the entangle- 
ment between spatial partitions of the same system, and also provides a distinct set 
of physical insights compared to the more standard spatial partitioning entanglement 
calculations. Particle entanglement provides a novel and unique perspective on the 
structure of itinerant many-particle wavef unctions. 
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Given a partitioning, the entanglement can be quantified using various measures. 
The basic quantity is the reduced density matrix of the A partition, pa = trg p, obtained 
by tracing out B degrees of freedom. We assume the system to be in a pure state, 
described by density matrix p = \ip)(%p\. Various entanglement quantifiers can then be 
extracted from pA- We will mostly confine ourselves to the entanglement entropy Sa, 
defined as Sa = — tr [p^lnp^]. We will also restrict ourselves to zero temperatures, i.e., 
to entanglement in the ground state of itinerant systems. 

Brief history - Several pieces of work explored simple versions of particle 
partitioning entanglement between identical quantum particles, even before the concept 
was carefully distinguished from spatial entanglement [1-7]. The relationship between 
quantum indistinguishability and entanglement was studied for two fermions in Ref. [1] 
and for two bosons in Ref. [2]. Refs. [4,6] studied particle entanglement in somewhat 
more complicated systems. 

A careful distinction with spatial entanglement, and a comparison between the two 
types of partitioning, appears in Refs. [8,9], in the context of fractional quantum Hall 
(FQH) states. Particle partitioning is tempting in entanglement considerations for FQH 
states, because FQH model wavefunctions (e.g., Laughlin states) are often written in 
first-quantized form where the particles have explicit labels. Thus, Refs. [4, 7] have also 
computed particle-partitioning entanglement in FQH states. 

In Refs. [8,9] and in work reported since then [10-12], particle entanglement has 
been shown to be a promising novel measure of correlations. In fractional quantum Hall 
states this type of entanglement reveals the exclusion statistics inherent in excitations 
of such states [8,9]. Similar insight arises from particle entanglement calculations in the 
Calogero-Sutherland model [11]. For one-dimensional anyon states, particle-partitioning 
entanglement is found to be sensitive to the anyon statistics parameter [12,13]. 

This review — Clearly, entanglement between particles in itinerant systems is a 
promising new concept, potentially useful for describing subtle correlations and the 
interplay between statistics and interaction effects. A broad study of the concept and 
its utility is obviously necessary. In this review, we will survey the results that are 
available until now. We will focus in particular on common features and on results of 
wide generality, that provide insights into classes of quantum itinerant systems. 

The present review is solely concerned with particle-partitioning entanglement in 
itinerant many-particle systems. We will therefore not discuss entanglement between 
spatial partitions, or any other kind of entanglement. Reviews of other types of 
entanglement can be found, e.g., in Refs. [14,15]. The other reviews of this special issue 
provide more recent and more condensed-matter-oriented perspectives on entanglement 
in many-particle states. 

The target audience for this review is condensed matter physicists interested in 
various possible kinds of entanglement in many-particle wavefunctions. As such, other 
than using the definition of the entanglement entropy Sa, we do not treat or use any 
quantum information theory topics. Ref. [16] provides a recent review of quantum 
entanglement from that perspective. It is of course not possible to make a review of the 
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present type completely self-contained, since we cannot introduce in detail each of the 
several models and states considered here. We therefore assume familiarity with several 
classes of many-particle states or models. Only minimal motivational background is 
provided for each model. The topic here is particle-partitioning entanglement and not 
the individual models. We expect that the typical practicing condensed matter theorist 
will indeed be familiar with most of the many-particle models and states employed. 

We start in Section 2 by working out in detail an elementary example of particle 
partitioning, contrasted to spatial partitioning of the same quantum state. The fact that 
entanglement depends crucially on the type of partitioning is perhaps not as widely 
appreciated as it should be; we hope a detailed example helps clarify the concept of 
particle partitioning. In Section 3 we present some generic results and intuitions, before 
moving on to specific systems in the following sections. 

Sections 4, 5, 6 consider respectively bosons, fermions and anyons, and review 
numerical and analytical results in order to provide an overview of various mechanisms 
for particle-partitioning entanglement. We then turn to more unusual many-particle 
states: section 7 reviews results for fractional quantum Hall states and section 8 for the 
Calogero- Sutherland models. 



2. An elementary example 



The concept of particle partitioning causes enough confusion to justify using a very 
simple example to illustrate in detail the definition and its difference from spatial 
partitioning. Readers comfortable with the concept may safely skip this section. 

We imagine two (spinless) fermions in three sites (or orbitals), which we label a, 
(3, 7. We will use the wavefunction 

\ip) = [aic^p + G^cjjcjyJ (vacuum) = a\ 1 1 1 ) + a 2 1 1 1 ) . 

For usual spatial partitioning, we can consider for example partition A to consist of site 
a only. Then the reduced basis for A consists of the a microstates |0) and |1 ) = cJjO), 

N 2 

|aip 

One could also take A to contain sites a and /3. Then the reduced basis for A 
contains four states, 1 00 ) , |01) = ct |0), |10) = c£, |0), and |11) = c^ct |0), and in this 



and the reduced density matrix in this basis is pa 



basis 
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/0 







\a 2 \ 





4|0), |10) 
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M 2 J 



We next turn to particle partitioning, for which the wavefunction must be expressed 
in first-quantized form with explicit anti-symmetrization: 

= a x [MVM 2 ) ~ </>a{2)M^)] + a 2 [^(1)0,(2) -^(2)^(1)]/^. 
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The particles (fermions) now have labels, so that we can consider the entanglement 
between particle 1 and particle 2. (Partition A contains particle 1.) Since particle 1 
can be in any one of the three sites, the reduced basis for A can be labeled by the site 
labels, \a), \j3), The reduced density matrix is 



From this simple example we already see that particle partitioning entanglement is 
utterly different from spatial or site partitioning entanglement. A second lesson 
that emerges from this example is that particle partitioning entanglement is affected 
substantially by the (anti-)symmetrization which is explicit in first quantization. It is 
therefore no surprise that this type of entanglement is especially sensitive to quantum 
statistics. 

3. General considerations 

Before analyzing specific systems, we present some facts and conjectures broadly 
applicable to a variety of itinerant quantum many-particle states. 

3.1. Bounds 

A generic itinerant lattice system has N particles in L sites; we consider bosons or 
spinless fermions with N < L. In every case, a natural upper bound for S n is provided 
by the (logarithm of the) size of the reduced density matrix pa = p n , i-e., the dimensions 

of the reduced Hilbert space of the A partition. This size is | ) = C(L,n) for 

V 71 ) 

fermions and C(L — 1 + n, n) for bosons. The actual rank of p n can be much smaller due 
to physical reasons, so that the entanglement entropies are usually significantly smaller 
than the upper bounds, as we shall see in the examples we treat. 

In a bosonic system, S n can vanish, since a Bose condensate wavefunction is simply 
a product state of individual boson wavef unctions, each identical. For fermions, however, 
anti-symmetrization requires the superposition of product states; for free fermions 
described by a Slater determinant wavefunction, this causes p n to have C (N, n) equal 
eigenvalues. This provides a nonzero lower bound for S n in a fermionic system. To 
summarize: 



PA 



( |ai| 2 /2 -aia* \ 
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Bosons : 
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3.2. Standard form for fermions 

For large fermion number, N^>1, we propose the following widely applicable form for 
the entanglement of n-CiV fermions with the rest: 

S n (N) = \nC(N,n) + a n + 0(1/ N^ 1 ) (3) 
= nlnN + a' n + 0(l/iV 7 ) , (4) 

with 7 > 0. This form is suggested by results reported in Refs. [8-12]. For example, 
a n = nlnm for the Laughlin state at filling v = 1/m [8]. The same standard 
behavior seems to hold for bosonic systems which lack macroscopic condensation into 
a single mode, e.g., bosonic Laughlin states [9], or hard-core repulsive bosons in one 
dimension [21]. Note that, for lattice sizes larger than N, the generic behavior (3) 
indicates that the entanglement entropy does not saturate the upper bound (1) or (2) 
obtained from the size of the reduced Hilbert space. 

Subtle correlation and statistics effects can be contained in the behavior of the 0(1) 
term a n , and sometimes also the Oil/ TV 7 ) term. To get some intuition about how such 
effects show up in a n , as we summarize the behavior of a n for several kinds of states. 
For free fermions, for charge-ordered states of the spinless-fermion model (subsection 
5.1, also Ref. [10]), and for Laughlin states (section 7 and Refs. [8,9]), we have 

a n (FF) = 0, a n (CDW) = ln2, a n (Laughlin) = n In m . 

We note that states which are intuitively 'more nontrivially correlated' have stronger 
n-dependence in a n . This strongly suggests that the a n function is a measure of 
correlations in itinerant fermionic states. It is natural to conjecture that the linear- 
in-n behavior of a n is symptomatic of intricately correlated states like quantum Hall 
states, and that in generic itinerant states a n will have sub-linear dependencies on n. 

3.3. Exceptional case of macroscopic degeneracy 

One has to treat with care cases where the single-particle spectrum has a highly 
degenerate ground state with degeneracy larger than N. The case of fractional quantum 
Hall (FQH) states is one example that we will treat in detail in Section 7. For FQH 
states, the appropriate Landau level is immensely degenerate and only partially filled. 
Another example is a macroscopically degenerate "flat band" that is partially filled. 
Flat bands appear in the band structure of frustrated lattices (kagome, checkerboard, 
pyrochlore,...), where one of the energy bands can be dispersion-free and have the same 
energy for all momentum; hence the name 'flat'. (See, e.g., Ref. [17] for a discussion of 
itinerant systems in flat-band situations due to lattice frustration.) 

Denoting the degeneracy by + 1 (FQH notation), we have the upper bounds 

s^j^y 5„<i D (^ +i ) + "- 1 ), 

respectively for fermions and bosons, from Hilbert-space counting alone. This is very 
similar to the case of iV particles in L sites, above. (The lower bounds are the same 



Entanglement between particle partitions in itinerant many-particle states 



In C(L,n) 



A 



In C(N,n) 



-r lnC(A^+l,n) 

V 



In C(N,n) 



Figure 1. Lower and upper bounds for particle-partitioning entanglement in fcrmionic 
systems. Left: usual situation such as mobile fermions in a non-frustrated lattice. 
Right: fermions in a flat band or within a Landau level, with macroscopic degeneracy. 
Dashed arrows show the effect of turning on interactions. 



as before.) The difference is that for non- interacting particles in a flat band, S n can be 
equal to the upper bound, whereas in the usual case it is equal to the lower bound. In 
this case of macroscopic degeneracy, interactions reduce the particle entanglement from 
the upper bound, as opposed to the conventional situation where interactions raise the 
particle entanglement from the lower bound. This is illustrated schematically in Figure 
1. 



3.4- Generic mechanisms: Slater terms; fragmentation or absence of condensation; 
massive degeneracies 

We summarize here the mechanisms through which an itinerant many-particle system 
can possess particle partitioning entanglement. This provides intuition concerning what 
physical information is contained in the entanglement between particle partitions. 

An ideal bosonic system is fully condensed in a single mode, and therefore has 
S n = 0. Thus, for bosonic systems, particle entanglement is a measure of the deviation 
from Bose condensation. Interactions provide a simple mechanism for this - since 
an interacting Bose system is only partially condensed, it possesses nonzero particle 
entanglement. Another mechanism is low dimensionality; hard-core repulsive bosons 
living on a line (ID continuum; Lieb-Liniger model) do not condense. In addition, 
condensate fragmentation provides a second mechanism for nonzero S n . Fragmentation 
is not as exotic a phenomenon as commonly perceived; in fact it provides a unifying 
perspective to describe Mott phases of bosons in lattice geometries [18]. Finally, one can 
also have S n ^ through quantum-mechanical superposition of condensates in different 
modes, i.e., Schrodinger cat states. In Section 4, we will illustrate these situations 
through several examples. 

For fermionic systems, fermionic statistics already provides a contribution 
InC (N, n) to the particle entanglement entropy - this is the value of S n when the system 
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wavefunction is a single Slater determinant. The excess particle entanglement over this 
amount tells us how far one has to go beyond a single Slater determinant in order to 
describe the physics of the system. In other words, the excess particle entanglement 
reflects loosely the number of Slater determinant terms of similar amplitudes that need 
to be combined to produce the system wavefunction. 

The cases of macroscopic degeneracies require additional considerations. Particle 
entanglement can be large here due to the much larger Hilbert space available without 
paying an energy cost, so that ground state wavefunctions can involve much more than 
a single Slater determinant (fermions) or a few condensate modes (bosons). 

3. 5. Relationship with correlation functions 

The purpose of this subsection is to clarify the relationship between particle partitioning 
entanglement and more traditional condensed-matter quantities, namely, correlation 
functions. 

In the easiest case, n — 1, the one-particle entanglement entropy S\ can be obtained 
from the one-particle correlation functions, or the single-particle state occupancies. For 
example, for one-dimensional translationally invariant systems, S\ can be expressed 
through momentum occupation numbers. The momentum occupation numbers are 
Fourier transforms of the reduced density matrix: c(k) = L^ 1 dxpi(x) cos(2irkx/L), 
with L being the size of the system. Then 

S 1 = -J2 c ( k ) lnc ( k ) ■ 

k 

More generally, for n > 1 it is intuitively clear that the elements of the reduced 
density matrices are proportional to correlation functions. We write down the precise 
relationship for a continuum one-dimensional case: 

Pn{xn,y n ) = ^1 W f (3?i) ■ ■ .(j)\x n )(t>{y n ) . . . (t>{yi)\^) ■ 
where the 'vectors' x n , y n encode coordinates of n particles, and 

Pn(x n ,y n ) = J dZ n+ l...dZ N ^*(xi,...,X n ,Z n+ l,...,Z N ) 

x ^(yi,...,y n ,z n+1 ,...,z N ) . 
The modification to lattice cases or higher dimensions is obvious. 

4. Bosonic systems: role of condensation 

In this section we illustrate the interplay between particle partitioning entanglement 
and Bose condensation, through several example models. First, considering a two- 
site Bose-Hubbard model, we demonstrate nonzero particle-partitioning entanglement 
through condensate fragmentation and the formation of Schrddinger cat states. Next, 
consideration of a lattice boson model and the continuum Lieb-Liniger model reveal 
bosonic particle entanglement due to partial condensation and absence of condensation, 
respectively. 
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4-1. Toy model: two-site Bose-Hubbard 

Following Ref. [10], we consider N bosons on a two-site 'lattice', subject to a Bose- 
Hubbard model Hamiltonian. The Hamiltonian is 

H= -(b\b 2 + blk) + \u (SISlSiSi + blblhk) ■ (5) 

We label the iV-boson basis states by site occupancies, i.e., as \Ni,N 2 ) = \Ni, N — N±) . 

For U = 0, the system is a non-interacting Bose condensate, with each boson 
packed into the single-particle state ^ (|1) + |2)). In the U — >■ +00 case, the system is 
a Mott insulator, with half the particles in site 1 and the other half in site 2, \N/2, N/2). 
Such a state is simple in the 'site' basis (second-quantized wavefunction), but involves 
symmetrization in the 'particle' basis (first-quantized wavefunction), leading to nonzero 
particle entanglement entropy. 

Finally, the U — > —00 limit involves all particles in either site 1 or site 2. The ground 
state is a linear combination of these two possibilities, ^75 (|0, N) + \N, 0)), which for large 
N is a macroscopic 'Schrodinger cat' state. Such a state is somewhat artificial, because 
an infinitesimal energy imbalance between the two states will 'collapse' this state. For 
example, a 'symmetry-breaking' term of the form e b\bi, added to the Hamiltonian (5), 
would favor site 2 and destroy the cat state. The resulting state is a product state with 
zero particle entanglement. 

For the simplest case of two bosons, there is only one way of partitioning (n = 1), 
so the only S n is Si. We expect Si = at U = 0, and maximal entanglement Si = In 2 
for both 'Mott' state at U — +00 and the 'Schrodinger cat' state at U = —00. The 
Hilbert space is small; one can diagonalize the problem and calculate Si analytically as 
a function of U. One finds Si (17) = Si(—U), interpolating smoothly between zero and 
In 2 ~ 0.6931 in both positive and negative directions (Figure 2a). 

Figure 2a also demonstrates the fragility of the cat state by showing the effect of 
an eb\i>i term (dashed curve). There is no appreciable effect for U > 0, but for U < 
the cat state is destroyed and we get Si — > for U — > — 00. 

For N > 2 bosons, it is meaningful to study S n with n > 1. The n-particle reduced 
Hilbert space has dimension n + 1; the reduced-space basis states can be labeled by the 
number of A bosons in site 1. In the Mott state \N/2, N/2), only the diagonal elements 
of p n are nonzero and they are all equal; hence S n (U — > 00) = ln(n + 1). In the cat 
state, only two elements are nonzero, both on the diagonal; hence S n (U —> —00) = In 2, 
independent of n. Figure 2b demonstrates, via calculation from wavefunctions obtained 
by numerical diagonalization, that S n increases to ln(n + 1) and In 2 in the U — > ±00 
limits. On the negative side, the change occurs sharply (around U = — 2/7V); the ground 
state remains nearly unentangled between U = and U = —2/N. 

Both p n (U) and S n (U) can be understood in greater detail using available 
approximations for the two-site model [18]; a description is given in Ref. [10]. 

To summarize, in the two-site Bose-Hubbard model the Mott state for U > 
and Schrodinger cat state for U < both possess particle entanglement. The 




Figure 2. (a,b) Particle partitioning entanglement for bosons in the ground state of 
a two-site lattice model with on-site repulsion, (a) Two bosons. Solid curve is for the 
basic Bose-Hubbard model. Dashed curve illustrates fragility of 'cat' state via a eb\bi 
term (e = 0.1). (b) 1-particle and 2- particle entanglement entropies for N = 10 (1000) 
bosons in two sites. Note different units for positive and negative U. (c) Entropy of 
spatial entanglement between the two sites (N = 2). 



particle entanglement in the two regimes have different physical origins: the physics 
of fragmented condensation for U > and that of Schrodinger cat states for U < 0. 

Comparison with spatial entanglement - - In Figure 2c we plot the entropy of 
entanglement between the two sites, i.e., the 'spatial' entanglement. The spatial 
entanglement is zero in the Mott regime of U — > +oo, where the wavefunction is a 
product state in the second-quantized wavefunction. (Symmetrization plays no role.) 
In the Schrodinger cat regime of U — > — oo, the spatial entanglement is In 2 like the 
particle partitioning entanglement. In the Bose condensate regime of small U, the 
spatial entanglement between sites is large, and for N = 2 peaks at In 3 at some small 
negative U. Thus particle partitioning and spatial partitioning lead to very different 
entanglements, except for the U — > — oo limit. 



4-2. Hard-core bosons on one- dimensional lattice 

Having considered fragmentation and cat states through the two-site model, we now 
turn to imperfect or partial condensation. One way to access such a state is through 
the model of hard-core bosons on a ID lattice (forbidden multiple occupancy, U = oo) 
with attractive nearest-neighbor interaction V: 

H = - ^ (c\c i+ i + 4 +1 CjJ + V^nin i+l + U ^^(rij - 1) 

i i i 
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with U — > oo. We consider N bosons in L sites, subject to periodic boundary conditions. 
This is closely related to the spinless fermion model treated later in subsection 5.1. 

The point V = —2 has a 'simple' ground state known exactly [27]. This 
wavefunction is a symmetric equal-amplitude combination of all possible C(L,N) 
arrangements of bosons. The exact wavefunction can be exploited to yield [10] 

S n = i/nlniV + C(iV ) 

where v = N/L is the filling fraction. A natural interpretation is that the pre-factor 
represents the un-condensed fraction. Whether this is generic for bosonic systems with 
partial condensation remains an open question. 

Details for n = 1 — The one-particle reduced density matrix is diagonal in the 
momentum representation, and has values 

(Jfe|pi|Jfe) = (N — \)/[L(L — 1)] for k^O 
(k\p x \k) = {L-N+l)/L for k = 

In the limit L — > oo (with fixed filling v = N/L), the k = eigenvalue becomes 
macroscopic at the expense of the others, indicating off-diagonal long-range order [19] 
and partial condensation with condensate strength 1 — v. 

Other cases of imperfect condensation — Interacting bosons in three dimensions 
also have partial condensation; it would be interesting to calculate S n for such a model. 

4-3. Lieb-Liniger bosons 

The generic one- dimensional continuum boson model with repulsive contact interactions 
(Lieb-Liniger model [20]) does not possess Bose condensation, and instead the particles 
fill up individual-particle levels just as fermions do. (In the Bethe ansatz these levels 
are labeled by rapidities.) One can thus expect a leading S n ~ nlniV behavior as for 
fermions. Currently available evidence strongly suggests this to be the case. 

For n — 1, the results of Ref. [12] (reviewed in Section 6) allow us to infer a 
leading IniV for the large-interaction limit, also known as the To nks- Girardeau limit. 
In addition, unpublished numerical results indicate that the behavior S\ ~ IniV holds 
for any nonzero interaction [21]. 



5. Fermionic systems: anti-symmetrization and correlations 

Following Ref. [10], in 5.1 we use numerical calculations of the spinless fermion chain 
with nearest-neighbor interactions (t-V model), one of the basic models of correlated- 
electron physics, to illustrate particle entanglement in fermionic systems. Other systems 
are commented on in 5.2. 
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5.1. Spinless fermions in one dimension 

We consider N spinless fermions on an L-site chain with periodic boundary conditions: 



Through a Jordan- Wigner transformation, this model can be mapped to the anisotropic 
Heisenberg (XXZ) spin chain model with anisotropy parameter A = V/2. 

For V = (free fermions), the ground state is simple in terms of momentum- 
space modes: a Slater determinant of the N fermions occupying the N lowest-energy 
modes. The n-particle reduced density matrix has C(N,n) equal eigenvalues, so that 
S n = In [C(N,n)], independent of the lattice size L. 

Half-filling — For repulsive interactions at half filling (N = \V), this model has 
a quantum phase transition at V — 2, from a Luttinger-liquid phase at small V to a 
charge density wave (CDW) phase at large V. This mirrors the well-known transition 
between gapless XY and gapped Ising phases in the XXZ model at the Heisenberg 
point A = 1 [22]. 

For N = |L, the ground state and hence particle entanglement can be simply 
understood in the infinite-^ limit. The ground state is an equal superposition of 
two 'crystal' states, and each of them gives a separate contribution to the reduced 
density matrix. The reduced density matrix has rank 2C(N,n) and equal eigenvalues: 
S n = In [2C(N, n)]. In the notation of Eq. (3), the sub- leading term a n interpolates 
between a n = at V = and a n — > In 2 at V — > oo for half filling. The interpolation 
details depend on n and N. 

Figure 3 show S n (V) for half- filling, calculated from wavefunctions obtained by 
direct numerical diagonalization. The S n (V) function evolves from Sff = ln[C(iV, n)] 
to In [2C(N, n)} ~ Sff + 0.6931. For n > 1, there is interesting non-monotonic behavior 
in some cases. At present there is no simple picture of the non-monotonic behavior. 

We also see Schrodinger cat physics in the t-V model: the V = +oo ground state 
is a superposition of two CDW states of the form 1 101010...10) and |010101...01 ). The 
fragility of this cat state can be seen by adding a single-site potential, ecjci, or a 
staggered potential, e'^Cg^. The ground state then collapses to a single crystal 
wavefunction, and S n drops to In [C(N,n)] (Figure 3 top panel). 

Phase transition — The small-n particle entanglement entropies show no strong 
signature of the phase transition at V = 2, even after extrapolating to the N — > oo 
limit. This is not too surprising because the notion of space enters rather weakly in the 
definition of particle entanglement; thus S n is not too sensitive to diverging correlation 
length or large-scale fluctuations. It remains unclear whether sharper signatures appear 
for finite n/N (as opposed to n<CiV). 

Away from half-filling — For N^L/2, the behavior is qualitatively similar to the 
half-filled case, a n increasing from zero to an 0(1) value as V increases from zero to 
infinity (Figure 4a). However, there is no simple picture for the V — > oo limit. Also, 
a n (V) appears to be monotonic, perhaps because a n (V — > oo) is not constrained as in 
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Figure 3. n = 1, n = 2 and n = 3 entanglement entropy in half-filled t-V model 
(N = L/2). The free-fermion contribution In [C(N, n)] has been subtracted off. The 
n = 1 plot also displays the effect of a symmetry-breaking ec[ci term, with e = 0.1. 
Inset: position of the maximum as function of e. 



the half-filled (CDW) case. 

Note that, except for S n= i in the half-filled case, the particle entanglement never 
saturates the upper bound, In [C(L,n)], dictated by Hilbert space size. 

Negative V - - An attractive interaction causes the fermions to cluster. In the 
V — > — oo limit, the ground state is a superposition (cat state) of L terms, each a 
cluster of the N fermions. The cat state can be destroyed as in the positive-^ case. For 
half-filling with even N, the V — > — oo wavefunction yields Si = hiiV-|-hi2 (Figure 4b). 
There are 0(N~ l ) corrections for odd iV = L/2. 



5.2. Other fermionic systems 

Other than the one-dimensionless spinless fermion model (5.1) and the quantum Hall 
states (Section 7), we are not aware of further explicit calculations for fermionic many- 
particle systems. Ref. [6] calculates particle entanglement of the Hubbard dimer (2-site 
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Hubbard model), and the current authors have found that preliminary numerics on short 
Hubbard chains show behaviors analogous to the spinless-fermion chain. 

6. Hard-core anyons in one dimension 

In one and two dimensions, quantum indistinguishable particles need not transform 
under exchange as fermions or bosons: rather, a continuum of possible intermediate 
cases connects the boson and fermion cases. Particles with such intermediate statistics 
are known as anyons [23]. 

A Bethe ansatz solution is available for the anyonic many-particle continuum 
model with contact interactions, and has recently received renewed attention [24]. The 
Hamiltonian is 

N d 2 

i 1 i<j 

and the anyonic statistics imposes the condition 

--,Xi, x i+1 , ■■■) = exp [i(9 - Tv)e(x i+1 - ■ ■ , x i+1 ,x h ■ ■ ■) , 

on the many-body wavefunction. Here e(x) = (or 1) if x > (x < 0)and 9 is the 
anyonic parameter. For 9 = or 9 = it this Hamiltonian reduces to free fermions 
or Lieb-Liniger Bose gas correspondingly. The choice of periodic boundary conditions 
*&(xi, ■ ■ ■ , Xi + L, ■ ■ ■ , xn) = ■ ■ ■ , Xi, ■ ■ -Xn) constrains the anyonic parameter to 

be an integer multiple of 2ix/(N — 1). Here L is the system size. 

The exact solution of the model, has been exploited in Ref. [12] to calculate the 
entropy of particle-partitioning entanglement Si between n — 1 anyon and the rest, in 
the limit 7—7-00. In this limit it is possible to compute the one-particle momentum 
occupation numbers, c e N (j) — 1/L p±(x) cos(2-njx/L)dx. This in turns allows one to 
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Figure 5. (a) One-particle entanglement entropy, Sf, as function of N, for 9 = ir (dots 
fitted by solid line), 9 = ir/2 (squares fitted by dashed line), and for 9 = tt (diamonds 
fitted by dotted line). The fits are of form Si ~ In N + f(9) + k{9)/Vn. (b) The filled 
dots are f{9) obtained numerically. The dashed line is a sine fit. 
(Data: courtesy of Raoul Santachiara.) 



obtain the one-particle entanglement entropy: 

oo 

Sl{N)= - 40')ln4(j). 

j=-oo 

Here j indexes the momenta, kj = 2irj/L. For free fermions (9 = 0) there are A" equal 
non-zero momentum occupation numbers: c° N (j) = 1 /AT for — (iV— 1)/2 < j < (N—l)/2. 
This gives the expected result Sf =0 (N) = \nN. 

For nonzero 9 it is not possible to obtain a closed analytic expression for the 
c e N (j). Ref. [12] provides some asymptotic relations (j <C N) using the Toeplitz 
determinant form of the one-particle density matrix, but this is not enough to extract 
the entanglement entropy. We will therefore only review numerical results extracted 
directly from the exact Toeplitz determinant. A fit to this numerical data (Figure 5a) 
indicates that in the limit A^ ^> 1 entanglement entropy behaves as: 

Sf(N)^lnN + f(9) + ^. 

f(8) is an A^ independent function and describes the main dependence of the one-particle 
entanglement entropy on the anyonic parameter 9. Figure 5(a) shows some S\ values 
calculated numerically, and from this one can extract f{9) values, plotted in Figure 5(b). 
The extracted f{9) values fit extremely well to a sine function. Explaining this regular 
behavior of the sub- leading term f{9) remains an open and important problem. 

The fact that entropy scales as lnA^ with a prefactor 1 is expected at 9 = and 
9 = 7r from the arguments involving fermions and condensation-less bosons, discussed in 
previous sections. In this sense, the leading lnA^ behavior for intermediate values of 9 
is no surprise. However, from the point of view of the momentum occupation functions 
c 9 N (j), the prefactor 1 is not evident, as the c 9 N (j) functions are quite nontrivial. For 
bosons (9 — 7r), the leading In A" behavior is reported to persist for finite values of the 
interaction 7 [21]. One might therefore expect a leading In A" at all 9 and all 7 > 0, but 
this has not yet been seen through explicit calculation. 
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7. Fractional quantum Hall states 

7.1. Preliminaries 

The fractional quantum Hall (FQH) states have long fascinated the condensed-matter 
community due to their remarkable transport properties and the exotic nature of 
their quasiparticle excitations [28-31]. Recently there has been enhanced interest in 
FQH states with non-abelian statistics [29-31], due to the possibility of implementing 
quantum computation schemes topologically protected from decoherence [32]. 

The unusual features of FQH states have been notoriously difficult to characterize 
using traditional condensed-matter concepts such as local order parameters and n-point 
correlation functions. Therefore, using novel measures of quantum correlations, such as 
entanglement entropies inspired by quantum information theory, is an attractive idea for 
quantum Hall states. One aspect of fractional quantum Hall states, namely topological 
order [33], has been successfully probed using spatial-partitioning entanglement [8,9,34- 
36]. Here, we will describe particle-partitioning entanglement in FQH states, following 
mainly Refs. [8,9]. The most striking result is that particle entanglements are closely 
approximated by upper bounds whose expressions reflect the exclusion statistics inherent 
in FQH states. 

We will consider both abelian and non-abelian FQH states, focusing on a 
paradigmatic example of each class, namely, the Laughlin (L) states [28] and the 
Moore- Read (MR) (or pfaffian) states [29,37,38]. In planar geometry, the respective 
wavef unctions are given by 

=n^-^-) m e- E ° N2/4 

i<j 

*mr({*» = Pf ( —^-) 11^ - %-) m e- E<NV4 , 

with Pf denoting the antisymmetric Pfaffian symbol, and = Xi + it/i representing 
the coordinates of the i-th particle. The fermionic Laughlin states (odd m) describe 
the physics of the most prominent series of FQH states observed as Hall resistivity 
plateaus in transport measurements on two-dimensional electron gases in semiconductor 
heterostructures. The Moore-Read state with m = 2 is widely expected to describe 
the state causing an observed plateau at a Landau-level filling fraction with even 
denominator. Since quantum statistics plays an important role for particle-partitioning 
entanglement, we also consider bosonic FQH states. These have not yet been realized 
experimentally, but have been the focus of experimental proposals and efforts with 
laser-cooled trapped bosonic atoms. 

We will describe FQH states in a spherical geometry [42] . In this representation the 
fermions are placed on a sphere containing a magnetic monopole. The magnetic orbitals 
of the relevant Landau level are then represented as angular momentum orbitals; the 
total angular momentum is half the number of flux quanta, L = \N ( j > . The + 1 
orbitals are labeled either / = to or L z = —L to +L. For N particles at fractional 
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filling v = 1/m, one finds the interesting FQH states for = mN — S, where S is 
a finite-size shift. The Laughlin states appear at S — m while for the Moore- Read 
states S = m + 1. The "filling" acquires the usual meaning v = N/N^ only in the 
thermodynamic limit. The orbitals are each localized around a "circle of latitude" on 
the sphere, with the / = orbital localized near one "pole." 

7. 2. Summary of main results 

Appreciating that derivations involving FQH states are not readily accessible to readers 
unfamiliar with the quantum Hall literature, we summarize our main results in this 
subsection. The technical details are deferred to the remaining subsections. 

We consider the entropy of entanglement between ua particles of the state and the 
remaining ub = N — n& particles. (In this section, we prefer to display the subscript A 
explicitly, because there is a profusion of symbols to distinguish from.) 

For both Laughlin and Moore-Read series of states, one can consider how the 
structure of the FQH wavefunctions reduce the rank of the reduced density matrices 
p nA . Hence one can derive upper bounds 5'^ ound for the particle entanglement entropy 
S nA [8,9], which are reduced compared to the naive bound (S F ) obtained from the full 
reduced Hilbert space. For FQH states on the sphere, in the simpler cases {tia = 2) the 
rank reduction has a physical interpretation in terms of the SU (2) multiplet structure 
of the spectrum of p nA . 

For Ua <C N, these upper bounds in fact happen to be close approximations to the 
actual values. This is due to the fact that the nonzero eigenvalues of p nA are distributed 
relatively flatly (Figure 6). (The more flat the eigenvalue distribution is, the closer the 
entanglement entropy is to the upper bound InD dictated by the local Hilbert space 
dimension T>.) This is in sharp contrast to the exponential-like eigenvalue distributions 
well-known from spatial entanglement [9,39-41]. 

The upper bounds are logarithms of combinatorial quantities which reflect the 
exclusion statistics of quasiparticle excitations in FQH states [25]. In addition, these 
quantities distinguish between the physics of the Laughlin and the Moore-Read states 
- the fact that the leading correlations have a 2-body nature in the m = 3 Laughlin 
states and a 3-body nature in the m = 2 Moore- Read states, is reflected in the 1/N 
expansions of the approximations S'^ ound . 

7. 3. Reduced ranks and entanglement upper bounds; fermionic states 

For N fermions, ua particles in the A block, and the total number of orbitals given by 
N<j, + 1 = 2L + 1, the obvious upper limit S\ from Hilbert space counting is: 



In the FQH states the correlations are such that the particles avoid each other and the 
entropy is further reduced. To obtain a handle on this, one may reason as follows. The 




(6) 
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model FQH states in the Laughlin and Moore- Read series can be characterized as zero- 
energy eigenstates of a Hamiltonian penalizing pairs and/or triplets of particles coming 
to the same position. After tracing out the coordinates for the B set, the dependence 
on those in the A set is such that one still has a zero-energy eigenstate. However, the 
number of orbitals available to the A particles is larger than what is needed to make the 
model FQH state in the A sector, and one instead has a certain number of quasi-holes 
on top of the A set model state. The total ground state degeneracy for this situation 
has been studied in the literature: see Ref. [38] for the Laughlin and Moore- Read states 
and Ref. [43] for the Read-Rezayi [30] and non-abelian spin singlet (NASS) [31] series 
of non-abelian FQH states. 

Laughlin states — The iV-particle Laughlin state is realized on a total of + 1 
Landau orbitals, corresponding to = m(N — 1) flux quanta. The Laughlin state 
for ha particles would need N A = m(n A — 1) flux quanta; we thus have an excess flux 
of AN^ = — N A = m(N — ua)- This corresponds to the presence of n qn = AN^ 
quasi-holes over the ground state. According to Ref. [38] each of the quasi-holes has a 
number of ua + 1 effective orbitals to choose from, with bosonic counting rules (meaning 
that two or more quasi-holes can be in the same effective orbital). This gives a number 
of quasi-hole states equal to 

(n A + 1) + n qh - 1 \ I n A + n qh \ ( Nj, + l — (m— l)(n A - 1 
n qh ) \ nA ) \ nA 

leading to the following upper bound to the entropy Sa 

£bound = ln f + 1 - ( m ~ l )( n A ~ 1) \ ^ 

We remark that this expression has a clear interpretation in terms of exclusion statistics 
[25]: the counting factor in Eq. (7) gives the number of ways ua particles can be placed 
in + 1 orbitals, in such a way that a particle placed in a given orbital I excludes 
particles from orbitals V with \l — V\ < m. 

Moore-Read — For the fermionic Moore-Read states at v = 1/m, with m = 2, 4, . . ., 
we can reason in a similar way, with now = m(N — 1) — 1. As for the Laughlin 
states we have an excess flux of AJV^ = — N A = m(N — ua) but now the number 
of quasi-holes is twice this number due to the fact that the fundamental quasi-holes 
correspond to half a flux quantum. Thus, n q h = 2AA^. We now take from Ref. [38] the 
following result for the total quasi-hole degeneracy 

\p f n qh /2 \ / (n A -F)/2 + n qh \ ^ 

F=n A mod 2\ F )\ ^ / 

The logarithm of this expression gives us an upper bound S'^ ound as before. 

Expansion in N^ 1 — For the Laughlin states, for ua <§C iV we get from equation 
(7) for large N 

1 <rji — 1 

S F_ abound = -^^UAiuA-l) + 0(1/N") 
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We compare this to the m = 2 Moore-Read states: 

S* - S A °™ d = jp\ n A{n A - l){n A - 2) + . . . 

The leading deviation from S A is a 3-body term at order 1/N 2 . This result nicely 
illustrates the fact that the leading correlations in the m = 2 Moore-Read state have 
a 3-body character: the wave-function vanishes if at least three particles come to the 
same position. 

For m^2 the leading correlations do have a 2-body character, as for the Laughlin 
states: 

S F A - Sr nd = ±—n A {n A - 1) + . . • 
N m 

Other fermionic FQH sequences — The quasi-hole counting rules for the order-fc 
clustered spin-polarized (Read-Rezayi [30]) and spin-singlet (NASS [31]) states are all 
known in the literature [43] . They can be used to generalize the upper bounds on particle 
entanglement entropy given in this subsection to these more intricate non-abelian FQH 
states. 

Bosonic quantum Hall states - - We consider bosonic Laughlin states at filling 
fraction v = — with m = 2, 4, . . .. The naive upper bound to the entropy associated to 
placing n A bosons in + 1 orbitals is 



^ = ln 



N<p + n A 
n A 



The expression for S'^ ound remains unchanged from the fermionic Laughlin case, giving 
the following leading correction in a 1/N expansion 

S B _ ^bound = }_ nA{nA _ x) + (Q) 

For a bosonic Moore- Read state, with filling fraction v = l/m with m — 1, 3, . . ., 
the leading 1/N correction becomes 

1 ffi 1 

C.B _ abound = _ 1} + . . . (10 ) 

N m 

In the case m — 1 the leading correlations have 3-body character, leading to the 
vanishing of the leading 1/N correction. (Similar to the fermionic m = 2 MR states.) 

Multiplet picture for rank reduction — For n A = 2, we can get a simple picture of 
the reduction of the entanglement entropy (or of the rank of reduced density matrix) 
compared to the naive bound, through consideration of multiplet structures present in 
the eigenspectrum of p nA . For FQH states on a sphere, the n^-particle reduced density 
matrices p nA commute with the total angular momentum operators L^ A and L z nA of 
the selected n A particles. As a result, the eigenvalues of p nA are organized in SU{2) 
multiplet structures: an eigenvalue for total angular momentum L nA will be (2L„ A + 1)- 
fold degenerate. 
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Figure 6. Eigenvalues for the 2-particle reduced density matrix, plotted against their 
multiplicities, for one Laughlin and one Moore-Read state. 



For ua = 2 fermions, each having angular momentum L = \N ( j ) , the 2-particle 
states have total angular momenta Li — 2L — 1, 2L — 3, . . ., 1 (0), for L integer (half- 
integer), giving a total number of (2L + l)(2L)/2 states. This corresponds to the naive 
upper bound to the entanglement entropy: 

S nA =2 < In [(21 + l)(2L)/2] = In ( N<t> ^ A J . (11) 

Inspecting the explicit structure of the fermionic Laughlin states with m = 3,5,..., 
one finds that the eigenvalues corresponding to 2-particle states with L 2 = 2L — 1, 
2L — 3, . . ., 2L — (m — 2) all vanish. The reason is that the correlations in the 
Laughlin states are such that particles cannot come too close together. For example, if 
a first fermion occupies the I = orbital, localized near the north pole, the Laughlin 
wavefunction has zero amplitude for finding a second fermion in orbit als 1 = 1,1 = 2, 
...,/ = m — 1 . The highest possible value of the angular momentum of the two fermions 
combined is thus L 2 = L + (L — m) . The remaining number of non-zero eigenvalues is 
(2L + (2 — m))(2L + (1 — m))/2, leading to an improved bound on the entropy 5*^=2 

Sn A =2 < In [{2L + (2 - m)){2L + (1 - m))/2] (12) 

which is clearly the n& = 2 case of equation (7). 

For ha > 2, the multiplet structures are more complicated, and it is difficult to 
generalize the above argument, but one expects that the arguments of the previous 
subsection, relying on the quasihole degeneracy, are equivalent to the vanishing of one 
or more eigenvalue multiplets. 



7.4- Numerical results 

Reduced spectra — In deriving the upper bound 5'^ ound we relied on the fact that a 
certain number of eigenvalues of the reduced density matrix vanish. The bounds would 
be exact if all non-zero eigenvalues were equal, but there is some eigenvalue spread, the 
bounds overestimate the actual values for the entropies. 

Figure 6 plots the eigenvalues for the = 2 reduced density matrix for iV = 9 
particles on a sphere in the m = 3 Laughlin state (L = 12). The horizontal axis 
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Figure 7. Particle entanglement entropies. Laughlin (m = 3) state in upper panels, 
Moore-Read (m = 2) state in lower panels. Dots are numerical exact values, the dotted 
line represents S% (eq. 6) and the solid curves are the bound S^f (cq. 7 for upper 
panels and eq. 8 for lower panels). 



represents the degeneracy 2L 2 + 1 of the eigenvalues, in descending order. The eigenvalue 
at Z/2 = 2L — 1 = 23, with degeneracy 47, vanishes; the non-zero eigenvalues show some 
scatter around an asymptotic value. Due to this scatter the entropy S = 5.509 is 
somewhat lower than the upper bound S^ ound = 5.533. 

An important difference between the m = 3 Laughlin and the m = 2 Moore-Read 
states is the absence of vanishing eigenvalues for the 2-particle reduced density matrix. 
The eigenvalue distribution shown in Figure 6 (right) illustrates this point. 

In the m = 2 Moore-Read state, there are vanishing eigenvalues in the reduced 
density matrix of ua > 3 particles. The number of nonzero eigenvalues predicted by 
Eq. (8) agrees with numerical results. For example, for n& = 3 and N = 10 particles 
there are 770 nonvanishing eigenvalues, in agreement with Eq. (8). 

While the nonzero eigenvalues are not all equal, their distributions are quite flat, in 
sharp contrast to the near-exponential eigenvalue distributions for spatial partitioning. 

Entanglement entropy values — In Figure 7 we compare numerically computed 
particle entanglement entropies with the bounds described above. 

Note that for the m = 2 Moore-Read case, 5'^ ound coincides with for ha — 2 
(dashed and solid lines coincide in lower left panel of Figure 7), which can be inferred 
from equation (8) and is also reflected in the fact that none of the eigenvalue multiplets 
vanish in that particular case (Figure 6 right panel). 
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7. 5. Reduced density matrices and correlation functions 

Since the n^-particle reduced density matrices p nA are obtained by integrating out all 
but ua of the particles, one expects these matrices to be related to the n^-particle 
correlation functions. For n& = 2, Ref. [9] provides a detailed explanation of this 
relationship, more specific than the general discussion of subsection 3.5. 

In particular, the n& = 2 eigenvalue distributions are indeed very closely related to 
the two-particle correlation function <?2(0- m fact, the g2{r) curve may be regarded as 
simply a continuous form of the A/ versus descending- (2/ + 1) curves of Figure 6. This 
explains why the eigenvalue distributions of p nA =2 in Figures 6, although discrete, are 
reminiscent of the well-known g2{r) curves for Laughlin and Moore- Read states. 



8. Calogero-Sutherland model 

The Calogero-Sutherland model holds a special place in condensed matter theory as 
an exactly solvable model which possesses fractional excitations [44]. The model 
has Jastrow-type ground state wavefunctions similar to FQH wavef unctions; the 
eigenfunctions are known in terms of Jack polynomials. Using properties of Jack 
polynomials ('duality relations'), Ref. [11] has studied particle partitioning entanglement 
entropies for this model. 

The model is described by the Hamiltonian (0 < x < L): 

A = jfl , y /^-l)(f) 2 

2 tt dx " ^ Sin2 (f & - ' 

By utilizing properties of Jack polynomials in the limit iV — n — > oo, Ref. [11] provides 
the following upper bound on the n-particle entanglement entropy S n in the iV-particle 
Calogero-Sutherland ground state: 

P(N-n) +n 
n 



S n < S b n ound = In 



In close analogy to the FQH results discussed in Section 7, this expression can be 
interpreted in terms of exclusion statistics. Tracing out one particle is equivalent to 
removing one particle from Fermi sea, which in turns creates /3 quasiholes. Thus taking 
out N — n particles creates (3(N — n) quasiholes. There are C(/3(N — n) + n, n) ways 
to accommodate n particles and /3(N — n) quasiholes within the Fermi sea, which gives 
the above estimate. 

Ref. [11] also gives an explicit expression for the sub-leading correction to this 
bound: 



S N ~ In 



In this formula the first term originates from the number of nonzero eigenvalues 
whereas the second one comes from the asymptotic eigenvalue distribution. An explicit 
expression for the O(N ) term is possible because of the more detailed understanding 
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of the eigenvalue distribution (in terms of the relevant Young tableaus) that is available 
for this model, as compared to the FQH case. 



9. Conclusions 

Particle entanglement is an emerging important measure of correlations in itinerant 
many-particle quantum systems. In this review, we have provided an extremely simple 
explicit example (section 2) designed to clarify the concept and remove misconceptions. 
We have synthesized the available results into a set of general observations (section 3). 
We have pointed out several different mechanisms for particle entanglement in itinerant 
quantum states, such as absent or imperfect Bose condensation, anti-symmetrization of 
fermionic systems, and Schrodinger cat states. 

The bulk of this review (section 4 onwards), of course, surveys the available results 
on particle entanglement in itinerant many-body systems. The discussion of more 
conventional bosonic, fermionic and even anyonic systems, falls mostly within the general 
observations of Section 3, with more subtle effects (e.g., anyonic statistics) showing up 
as sub-leading terms. The study of particle partitioning in quantum Hall states and 
the Calogero-Sutherland model provides a window into more exotic phenomena, such 
as exclusion statistics. 

Clearly, the study of particle partitioning entanglement is at its infancy, and one 
expects further insights and new calculations to appear. We hope this review will provide 
inspiration for further advances. 
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